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Fully-developed incompressible Navier-Stokes turbulence in three dimensions is a dissipative dy- 
namical system that exhibits strong departure from absolute equilibrium. Nevertheless, several 
kinds of representation by Tsallis equilibria have been proposed. The symmetry of the contribu- 
tions to the kinetic energy from the degrees of freedom of the flow must be broken in order to 
construct turbulence applications. Tsallis representations of turbulence involve the extrapolation 
of results for systems with a single degree of freedom, a procedure that requires critical examina- 
tion. Theories of Tsallis statistics for acceleration of fluid particles are compared with computer 
simulations of the pressure gradient field in the present paper. Applications of Tsallis formulations 
to statistics of longitudinal velocity differences are analyzed. No flavor of equilibrium statistics 
based on the kinetic energy can encompass the dynamically fundamental turbulent energy cascade, 
which requires correlation of triads of degrees of freedom at different scales. 
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1. INTRODUCTION 

The search for something that is maximized or mini- 
mized by turbulent flows extends back over more than a 
half-century. Non-extensive entropy has joined a varied 
set of candidates. 

Fully developed incompressible Navier-Stokes (NS) 
turbulence is a dissipative dynamical system with very 
many degrees of freedom in a state of strong departure 
from absolute statistical equilibrium. The disequilibrium 
is expressed by several features. In isotropic NS turbu- 
lence in three dimensions (3D), the wavenumber energy 
spectrum E(k) at high Reynolds numbers goes approx- 
imately like /c -5 / 3 in the inertial range of k, which ex- 
tends between the low k where any driving occurs and the 
high k where viscous dissipation is strong. In contrast, 
Boltzmann-Gibbs equilibrium based on kinetic-energy 
yields E(k) oc k 2 . Moreover, there is strong energy- 
transferring coupling within the inertial range between 
k that differ by a factor greater than 2. If one thinks of a 
temperature local in fc, this means strong coupling among 
degrees of freedom whose local temperatures differ by a 
factor like 2. 

In 1988 Tsallis proposed a generalization of statisti- 
cal mechanics that can describe non-extensive statistical 
equilibrium, in contrast to other generalizations, like that 
of Renyi, that preserve the additivity of the entropies of 
independent subsystems [1] . The Boltzmann-Gibbs equi- 
librium is recovered as a limiting case. Tsallis statistics 
have been formulated for a variety of systems including 
fluid turbulence. The theoretical formulations for 3D NS 
turbulence assume isotropy. The results have been com- 
pared not only with approximately isotropic direct nu- 



merical simulation (DNS) but also with experiments that 
have more complex geometries. The key analysis typi- 
cally is performed for systems very much simpler than 
turbulence: one-dimensional maps and dynamical sys- 
tems consisting of a single degree of freedom that obeys 
a modified Langevin equation. 



2. TSALLIS STATISTICS 

A number of generalized definitions of entropy have 
been proposed that reduce to Boltzmann-Gibbs entropy 
in a limit. The entropy proposal by Renyi [1] was moti- 
vated by information-theoretic considerations. It yields 
additivity of the entropies of independent subsystems. 
Tsallis statistics are based upon the entropy formula 
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where q is a parameter of the theory and the sum is over 
the probabilities pi of states i that assign values to a 
complete set of real dynamic variables y n . 

If q 7^ 1, S q is not the sum of contributions from each i. 
Independent subsystems give interlocked contributions to 
S q , yielding non-extensive equilibrium. The limit q — > 1 
yields the Boltzmann formula S q — —J2Pi^ n Pi- 

Let Ei be the energy of state i. The extreme value of 
S q under the generalized energy constraint 
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is realized by the probabilities 
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Pi = ^[l + (q-l)(3E i }- 1 ^- 1 \ 

Here k B (3 is an inverse temperature parameter whose 
value is determined by the value taken for U q , k B is Boltz- 
mann's constant and 

^ = X)[ 1 + («- 1 )^] _1/(,_1) - 

The Lagrange multiplier for the energy constraint is 
written in the form (q — 1)(3 so that the limit q — > 1 
yields the Boltzmann-Gibbs values 

Pl = exp(-PEi)/Z B , Z B =Y, eM-PEi). (5) 



Consider the case in which (a) each Ei is a sum of 
energy contributions e n from the y n ; (b) the dynam- 
ics conserve energy; (c) there is a Liouville property 
^2 n dy n /dy n = 0. In this case it is self-consistent to 
take uniform density of states in the phase space of the 
y n . The Boltzmann-Gibbs formula (5) then yields an ex- 
plicit result for the probability density function (PDF) of 

Vn- 



P n (y n ) = exp(-/?e n )/S„, S n 



cxp(-/3e„)dy n . 

(6) 



If q ^ 1, Ei that are sums of contributions from each y n 
do not yield separation of the PDFs of distinct y n . In 
most cases P n (y n ) can be obtained only by integration 
over all the y m ^n- 
The particular case 



Ei = [u 2 /2]i = 



n = 1, 2, ...N, 



is appropriate to a suitably constructed truncated rep- 
resentation of Eulcr or NS flow. In this case only, the 
constant energy surfaces are hyperspherical surfaces in 
the phase space of the y n . The Boltzmann formula (6) 
immediately yields Gaussian P n {y n ). 

If g ^ 1, (3) and (7) yield complete symmetry in n. 
An explicit result for P n can be extracted if, in addition, 
TV — > oo. Then the distribution of each y n on the energy 
shell tends to the Gaussian 

[Pn{y n )]i = -J^eM-lNyl/u 2 ). (8) 

U V ^7T 

With the use of (3), (8) can be integrated over u to give 



^ The canonical distribution (5) is realized if the con- 

servative system is put in contact with a heat bath and 
then removed. This implies unique stability properties 
for the canonical distribution that can be demonstrated 
as follows. Consider two systems A and B that are ini- 
tially independent with PDFs Pa(E a ) and P A (E B ) that 
(4) are functions only of the system energies. The initial 
PDF for the combined system is then Pa(Ea)Pb(E b ). 
Now let the systems be coupled by arbitrary couplings 
that conserve the total energy Ea + E B and preserve the 
Liouville property. The condition that the total system 
PDF be invariant under these couplings is that the PDF 
keeps the form Pa(Ea)P b (E b ) and is time invariant. 
But since there are arbitrary couplings, the PDF must 
also have the form Pa+b(Ea + E B ). Equality of these 
two forms has the canonical distribution as its unique so- 
lution, a fact easily verified by taking the derivative with 
respect to Ea and separately with respect to E B , then 
rearranging terms. 

Tsallis q ^ 1 statistics for a conservative system with 
Liouville property do not have the robust stability prop- 
erties of the canonical distribution. A necessary condi- 
tion for q 7^ 1 is that the system be isolated from its en- 
vironment or the environment be specially constructed. 

Dissipation destroys the Liouville property. Equilib- 
rium statistics can be extended to some dissipative sys- 
tems, such as a single degree of freedom that obeys a 
Langevin equation with suitable relation between damp- 
ing and random forcing. If particular assignments of ef- 
fective energy are made, dissipationlcss equilibrium Tsal- 
lis statistics yields exact PDFs for some dissipative sys- 
tems. Examples of such systems, with a single degree of 
freedom, will be discussed in Sec. 4. NS turbulence at 
high Reynolds numbers is a dissipative system with many 
strongly interacting degrees of freedom. In the forced 
steady state, the forced degrees of freedom are separated 
from those where the dissipation occurs by a multistage 
energy cascade. Any attempt to describe NS turbulence 
by equilibrium statistics must be considered courageous. 
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P n (y n ) 



f [1 + \P{q - iVr 1 ^ eM-^yl/u 2 )^. 
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3. APPLICATION TO NAVIER-STOKES FLOW 



An underlying fact is that the quadratic inviscid con- 
stants of motion give no hint of the rich dynamics that 
emanate from the NS equation. For example there is 
nothing to imply the crucially important Kelvin circula- 
tion theorem. 

The exploration of the Boltzmann-Gibbs and Tsallis 
statistics of Sec. 2 for NS flow begins with the choice of 
(9) a representation basis for the flow. There is an infinity of 
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representation bases, constrained by the boundary con- 
ditions on the flow. The DNS to be described in what 
follows uses cyclic boundary conditions on a cube. This 
is appropriate for spatially homogeneous flow. The rep- 
resentation most immediately suggested by homogeneity 
is expansion of the velocity field into wavevector compo- 
nents. 

In order that the system be finite, the expansion can be 
truncated at a wavenumber fc max sufficiently large com- 
pared to typical dissipation wavenumbers that excitation 
above k max can be neglected. The real variables y n at 
each k are the amplitudes of the real and imaginary parts 
of the flow components perpendicular to k. The kinetic 
energy is given by (7) if the fluid density is set to unity. 
Orthogonal transformations yield other representations 
in which Ei is a sum of squares. One class of such bases 
is wavelets, which may be appropriate to describe the 
intermittency of turbulent flow. 

High Reynolds number turbulence maintained in sta- 
tistically steady state by forcing at large spatial scales 
and viscous dissipation at much smaller spatial scales 
exhibits a cascade of energy through intermediate scales, 
the inertial range. Energy cascade is the most fundamen- 
tal dynamical feature of the turbulence. It is not a phe- 
nomenon of absolute equilibrium, whether Boltzmann- 
Gibbs or Tsallis. Energy cascade requires nonzero triple 
correlations of the form (y n ymyr), where n,m,r label 
modes with different wavenumbers or other scale labels. 
Such correlations are not produced by any flavor of equi- 
librium statistics based on the kinetic energy: the ki- 
netic energy is a sum of squares and is invariant under 
y n — > — y n . The characteristic times of the flow modes 
associated with energy cascade and other dynamical fea- 
tures decrease with scale size but are finite at all scales. 

In any representation where (7) holds, the results of 
Sec. 2 imply that neither Boltzmann-Gibbs nor q ^ 1 
Tsallis statistics describes high Reynolds number turbu- 
lence. Both statistics yield equipartition of kinetic energy 
among the degrees of freedom because of the symmetry 
of (7) in the y n on each energy shell. An application of 
Tsallis statistics to turbulence, viable in the sense of fit- 
ting certain modal PDFs, must break the symmetry over 
the degrees of freedom. In some work, individualized val- 
ues of q and (3 are invoked for subsets of the degrees of 
freedom and/or the kinetic energy is replaced by a modi- 
fied energy function that is different for different degrees 
of freedom. 

Different q and values for two subsets of the flow 
degrees of freedom may possibly be justified as an ap- 
proximation if the subsets have dynamical coupling suf- 
ficiently weak compared to the strength of their internal 
dynamics. If the subsets are ranges of wavevector am- 
plitudes in a homogeneous flow, both approximate the- 
ory and DNS suggest that this requires separation of the 
ranges by more than a decade in k. If some weakness of 
coupling criterion is not satisfied, there is no justification 



for assigning different q and j3 values to different degrees 
of freedom. 

The turbulence PDFs that have been approximated 
in applications of generalized equilibrium statistics have 
shapes that are easy to fit in a number of ways: the 
PDFs are smooth and nearly symmetrical with flaring 
skirts. Justification of a fit from turbulence dynamics is 
essential for the fitting process is to be considered mean- 
ingful. 

The applications by Arimitsu and Arimitsu use the link 
between Tsallis statistics and multifractals [2,3] to deter- 
mine a single q < 1 value for the entire system, related to 
the intermittency behavior associated with a multifrac- 
tal turbulence model. The multifractal model describes 
intermittency that increases with decrease of scale size. 
Multifractal models of turbulence have not been derived 
from the NS equation but they are supported by theo- 
retical arguments and their parameters can be tuned to 
agree well with a variety of experimental measurements. 

Many of the applications deal with velocity increments 
across spatial separations. Velocity increments do not 
constitute a representation basis for a flow. In contrast 
to the invertible transformation between, say, wavevector 
components and wavelet amplitudes, the transformation 
from the velocity field in x space to sums and differences 
of velocities at two space points is overcomplete. If the 
velocity sums are struck out, leaving only the increments, 
the transformation is not invertible. Analysis that de- 
rives results by working with velocity increments rests on 
the implication of corroboration in a representation basis 
from which the velocity increments are constructed. 

In Sections to follow, we shall examine several appli- 
cations of Tsallis statistics to turbulence. Before that, 
we shall summarize the behaviors of three systems with 
a single dynamical variable for which Tsallis statistics is 
exactly realized: anomalous diffusion in an optical lattice 
[4] , a Langevin equation with fluctuating parameters [5] 
and sensitive nonlinear maps [2,3]. 



4. TSALLIS STATISTICS FOR SINGLE 
VARIABLES 

Lutz [4] finds that the diffusion of momentum p in 
an optical lattice treated semiclassically satisfies the 
Langevin equation 

dp/dt = K(p) + dD(p)/dp + y / 2D{p)w{t) (10) 

where K(jp), D(p), and dD{p)/dp are nonlinear non- 
stochastic functions of p and w(t) is Gaussian noise with 
(w(t)w(t')} = 5(t - t'). The exact PDF in steady state 
that follows from (10) is 

P p (p) = Z- 1 [l + (3(q-l)p 2 }- 1/{q - 1) . (11) 
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Here q > 1 and (3 are found from D(p) and K(p); p 2 
comes from the specific forms of K{p) and D(p); Z is the 
normalization factor. 

Lutz notes that the steady state described by (10) and 
(11) is not an equilibrium state and (3 is not a well-defined 
inverse temperature. 

Beck [5] has demonstrated that Tsallis statistics with 
q > 1 are realized by a dynamical system consisting of a 
degree of freedom that obeys a Langevin equation with 
fluctuating parameters. Let y satisfy 

dy/dt + ly = <jw(t), (12) 

where 7 is a damping constant, a is a strength parame- 
ter, and w(t) is unit white-in-time Gaussian noise. The 
steady-state PDF of y is Gaussian with variance 1/(3, 
where (3 — j/cr 2 - 

Now suppose that 7 and/or a fluctuate so that (3 = 
7/cr 2 has a x 2 distribution of degree n. This means that 
the PDF of (3 is oc /3"/ 2 " 1 exp(-n/3/2/3 ), where O is the 
mean value of {3. Assume that (3 fluctuates on time scales 
infinitely large compared to the 0(1/7) time required for 
y to reach equilibrium. Beck shows that the PDF of y 
has the form 

P{y)^l/[l + P'(q-l)y 2 ]- 1 '^- 1 \ (13) 

where 

g=l + 2/(n+l)>l, (3 1 = 2/? /(3 - q). (14) 

In comparison with the formulas of Sec. 2, y 2 plays the 
role of an energy. 

Although (10) and (12) are dissipative equations, the 
PDFs (11) and (13) express Tsallis non-dissipative equi- 
librium if the energy constraint is constructed with the 
required effective energies. 

Certain classes of nonlinear one-dimensional maps have 
a kind of sensitivity to initial conditions, at the onset of 
chaos, that can be described by Tsallis statistics. The 
growth of small perturbations in initial conditions is not 
exponential but instead follows a power law that corre- 
sponds to a Tsallis distribution with 1 > q > 0. Lyra 
and Tsallis [2] find that the critical attractor for several 
kinds of maps, including generalized logistic maps and 
the circle map, exhibits multifractal behavior in which 
the singularity strength a has maximum and minimum 
values that are related to q in a Tsallis distribution by 

1/(1 - q) = l/a min - 1/ttmax- (15) 



5. ACCELERATION OF FLUID ELEMENTS 

The statistics of pressure forces, viscous forces, and the 
resulting acceleration of fluid elements in incompressible 



NS flow has drawn continuing attention for more than 50 
years [6-9]. The behavior of fluid-element acceleration 
(acceleration along Lagrangian trajectories) can be ob- 
tained completely if the distribution of the velocity field 
is taken to be Gaussian with a specified spectrum. Sup- 
pose that the Gaussian velocity field is isotropic, with a 
wavenumber spectrum that includes a long fc -5 / 3 range 
followed by a rapidly falling dissipation range. Then 
the acceleration a is dominated by pressure forces, which 
make a contribution 

(a 2 ) = c^/ 2 ^ 1 ' 2 (16) 

to the mean-square (a 2 ), where c p is a number nominally 
0(1), e is the average dissipation per unit mass, and v is 
the kinematic viscosity. The viscous forces make a con- 
tribution of the same form and nominal magnitude, but 
the corresponding coefficient c v is smaller than c p by a 
large enough factor that a is strongly dominated by the 
pressure-gradient contributions. 

In realizations of approximately isotropic NS turbu- 
lence at large Taylor microscale Reynolds number R\ by 
DNS or experiment, intcrmittency strongly affects the 
distribution of pressure-gradient forces and causes a small 
deviation from the —5/3 energy spectrum exponent for 
the inertial range. The effect on the pressure-gradient 
contribution to (a 2 ) probably is to make c p weakly depen- 
dent on the ratio kd/ko of the dissipation and macro-scale 
wavenumbers. It is to be expected that the dominance of 
pressure over viscous contributions to a is enhanced com- 
pared to the Gaussian case, especially for contributions 
to the tails of the PDF of a. 

Beck [10] constructs Tsallis statistics for the acceler- 
ation, with the conclusion that the normalized PDF of 
x = a/(a 2 ) 1 / 2 is 



The analysis uses a single variable u defined as the dif- 
ference in the velocity measured at two times along a 
fluid-element trajectory. The time difference r is taken 
to be the order of the characteristic dissipation time 
t v = (v/e) 1 ' 2 of the 1941 Kolmogorov theory (K41). The 
acceleration is measured as o = u/t. 

Beck models the dynamics of u by that of y in the 
Langevin equation with Gaussian and \ 2 fluctuating pa- 
rameters whose analysis is summarized in Sec. 4. He 
states that an excellent fit to experimental data on the 
PDF of a [11] and on (a 2 ) [12] is obtained by taking 
q = 3/2 and using, instead of (3' in (13), an O^" 1 / 2 ^" 1 / 2 ) 
parameter whose precise value is tuned. The choice 
q = 3/2 gives (17). 

In more recent work (private communication and [10]), 
a fit to normalized pressure gradient component data 
x i = Pi/ipf) from approximately isotropic DNS at R\ = 
381 [13] is made using log-normal rather than x 2 statis- 
tics in (12). There is one fitting parameter. The bold 
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simplicity of (17) is lost, but the fit is extended over the 
range to P x (x) ~ 10~ 12 . The \ 2 an d log-normal fits 
are plotted with the DNS PDFs in Figs. 1 and 2. Fig. 
2 shows that neither fit is accurate in the central peak 
region, which contains most of the weight of P x (x). 

A fluctuating equilibrium environment, of which the 
log-normal and \ 2 statistics are particular examples, is 
inconsistent with the turbulent energy cascade. See Sec. 
8 for further discussion. The white-noise forcing in (12) 



conflicts with the fact that turbulence modes of given 
scale are principally driven by modes of larger scale that 
have finite and longer characteristic times. 

The DNS pressure gradient curves agree well with ex- 
perimental acceleration data [11] measured in the approx- 
imately isotropic central region of flow between counter- 
rotating disks at reported R\ values of 200, 690, and 970. 
The experimental data extend down to about 10~ 6 on the 
vertical axis. 
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FIG. 1. Lin-log plot of normalized PDFs of components of DNS pressure gradient at R\ — 381 together with Beck's \ 2 an( l 
lognormal fits. The tails of the DNS PDF curves fit a stretched exponential of form exp(— |a;| 0,41 ). 



Arimitsu and Arimitsu [3] derive an acceleration PDF 
from a multifractal cascade model. The starting point 
is (15), which they generalize to the form 1/(1 — q) = 
\joL- — l/a + , where negative f(a) is admitted and a_ 
(a+) is the smaller (larger) zero of the multifractal spec- 
trum f(a). The spectrum /(a) is linked to that of a 
cascade model for isotropic turbulence. Specifically, 

f(a) = 1 + (1 - q)- 1 log 2 [l - (a - a ) 2 /(Aa) 2 }, (18) 

where (Aa) 2 = 2X/[(1 — q) In 2] and q, X, ao are deter- 
mined from the intermittency exponent /i. The multi- 
fractal cascade is assumed to proceed in decrement steps 
of 1/2 in "eddy size". After m statistically independent 
steps the probability P^(a)da to find an eddy with 
singularity in the range a to a + da is determined as 
The acceleration statistics arc obtained from the scal- 



ing 6p m /5p = {l m jk 



\2a/3 



Here Sp m is the pressure 



difference across the distance l m = 2~ m £ and £ is a tur- 
bulence macroscale. The cascade is terminated at step 
n chosen so that Sp n /£ n is expected to approximate the 
pressure gradient. 

An additional element in the analysis is "thermal fluc- 
tuations" that have a Gaussian or q > 1 Tsallis distri- 
bution and are expected to dominate the acceleration 
PDF for small amplitudes. The total PDF is written as 
the sum of a contribution from the cascade up to step 
n and the "thermal fluctuations" part. We do not find 
a physical basis for the "thermal fluctuations" . If such 
fluctuations existed, we would expect them to pervade 
space so that they add everywhere to the field produced 
by the multifractal cascade. In that case, the PDF of the 
total field would not be the sum of PDFs from the two 
contributions to the field. 
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FIG. 2. Same as Fig. 1 but only the central regions of the PDFs are shown on a lin-lin plot. 
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FIG. 3. Lin-log plot of PDF of normalized DNS pressure gradient component Px/iPx) 1 ^ 2 a ^ R\ ~ 381 compared with the 
PDF proposed for fluid-element acceleration in [3]. Curves for cascade termination step numbers of n — 20.9, 9.0, and 5.0 are 
shown. For these fits, fi = 0.210 and q = 0.3183. 



P {m \a) = [pW] ro a [1 - (a - a ) 2 /(Aa) 2 ] m/(1 ~ 9) - [ 



(19) 

We have calculated acceleration PDFs from the com- 
bined multifractal and Gaussian theory as presented in 



3]. The results are compared with DNS [13] of pres- 
sure gradient in Figs. 3 and 4 for n — 20.9, n = 9, and 
n = 5. The curves for these n values are obtained with 
H = 0.210, q = 0.3183. R\ = 381 in the DNS. 

The DNS is cyclic in each axis direction and the num- 
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ber of points in each direction is 2 10 = 1024. The largest 
possible value of £q pointing along an axis therefore is 
2 9 = 512 grid point spacings (this is a distance n). Thus 
the upper limit for n is 9. In fact, this limit should 



be reduced because the forcing extends from the lowest 
wavenumber k = 1 to k = ^/6. The multifractal cascade 
consequently should be taken to start at no more than 
512/V6 w 209 grid point spacings. 



Oh 



0.8 



0.6 



0.4 



0.2 







Gaussian 
n=20.9 
n=9.0 
n=5.0 
DNSp x 




-3 



■1 



x 



FIG. 4. Same as Fig. 3 but only the central regions of the PDFs are shown on a lin-lin plot. 



There is substantial scale nonlocality in the turbulence 
coupling so that there must be a non-asymptotic region 
at the start of the assumed multifractal scaling range. If 
we assume that this region starts at approximately 209 
grid point spacings and extends for approximately 3 of 
the nominal cascade halving steps, we get 5 as a plausi- 
ble value of n. If R\ were increased to 1000, the inertial 
range would be extended by approximately two scale-size 
halvings. 

Although the value n = 20.9 is not meaningful for the 
DNS, we show it because it gives optimum fitting of the 
theoretical curve to the skirts of the DNS curve. We find 
that the optimum fit obtained by tuning both [i and n 
is insensitive to variation of fi over the range 0.210 to 
0.260 and a corresponding variation of n from 20.9 to 
17.0. Over this range the resulting theoretical curve is 
almost invariant. There appears to be no way to bring 
the tails of the theoretical curve for n = 9.0 or n = 5.0 
close to the DNS by tuning /j,. 

The best fit stated in [3] uses fj, = 0.241, n = 18.2. 
This is almost identical with one of the fits we found 
(/i = 0.240, n — 18.0) in the parameter domain where 
the result is nearly invariant to fitting changes. Note 
that n = 18 would correspond to 1/512 grid point spac- 
ing if £o is given the maximum possible value of 512 grid 



point spacings. 

The crucial parameter in the fits is the number of sta- 
tistically independent cascade steps, which is the maxi- 
mum value n of m in (19). In the preceding discussion, 
we have started with a maximum meaningful value of £ 
and assumed that the scale decrement in each cascade 
step is 1/2. Then n = 18 yields the meaningless min- 
imum scale size 1/512 grid point spacing. Arimitsu & 
Arimitsu discuss two alternatives to force the minimum 
scale size to be the order of one grid point spacing with 
n = 18. One is to make £ much larger than the largest 
meaningful value of 512 grid point spacings. The other 
is to make the inverse of the scale decrement factor less 
than 2. It must be stressed that a factor of 2 is already 
is smaller than the scale ratio needed for approximate 
dynamical independence of flow modes. 



6. SCALE DEPENDENT TSALLIS PARAMETERS 

K41 implies that velocity statistics associated with the 
succession of scales in the inertial range are similar. It 
has long been known from experiment and suggested by 
theoretical arguments that this is not the case. Instead, 
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there is intermittency that increases with decrease of spa- 
tial scale. One manifestation of the intermittency is the 
behavior of the velocity structure functions, defined as 
moments of order n of velocity differences across a dis- 
tance r. As r takes successively smaller values in the in- 
ertial range, the structure functions show an increasingly 
rapid growth with n, in conflict with K41 similarity. The 
growth of the structure functions expresses anomalous 
scaling exponents. A number of theoretical models have 
been explored, but almost nothing has actually been de- 
duced from the NS equation. 

Beck has studied a series of models of the longitudinal 
velocity difference u (difference in the velocity component 
pointing along r). We shall consider them in more or less 
chronological order. In the first of these models [14], the 
dynamics of u is related to a Langevin equation with 
non-stochastic damping and a stochastic forcing term: 



du/dt + ju = f(t). 



(20) 



Here f(t) is what Beck calls deterministic chaotic noise 
with typical time scale r <C 1/7. The statistics of / is 
such that the PDF of u acquires a skewness. The actual 
qualitative physics of transfer into modes of given scale 
is that the transfer comes mostly from modes of larger 
scale with longer dynamical times. This conflicts with 
r <C I/7 in the model equation (20). 

The result is that for c/y/E^ < l/u the PDF of 
u = uj \J (u 2 ) has the form 



l + 0q 



¥ 2 



(21) 



where q is a Tsallis entropy parameter, q = q — 1, R\ 
is the Taylor microscale Reynolds number of the turbu- 
lence, and c is a positive O(l) constant. Higher-order 
terms in c prevent a change of sign of p(u) at large nega- 
tive u and make p(u) normalizable. The value of (3 is ad- 
justed to give u variance 1, which requires (3 = 2/(5 — 3q) 
for large R\. 

How is (21) obtained from (20) and how is q deter- 
mined? There are a number of steps. The small pa- 
rameter jt is identified with 1/R\ by the following ar- 
gument: since 7 is a damping, it should be proportional 
to the underlying viscous damping coefficient v\ correct 
dimensionality then requires a length scale, which is set 
equal to the Taylor microscale A, so that 7 = v/\ 2 . Also 
it is assumed reasonable to take r = A/ \J (it 2 ). Thus 
^/yr = c/Vi?A with c = 0(1). 

The parameter q — 1 is set according to the number of 
cascade steps between r and the K41 dissipation length 
7] = (^ 3 /e) 1 / 4 . A cascade step is taken to be a factor of 
2 in 1/r. This leads to 



l/(g-l) = 0(l)+Iog 2 (r/ >7 ). 



(22) 



Finally O(l) is set to 1 because that is said to give agree- 
ment with experimental data. 

Beck uses (21) to calculate moments of the normalized 
velocity difference u (normalized structure functions) and 
from them calculates scaling exponents for the normal- 
ized structure functions. 

Some comments on (20)-(22): In K41, the characteris- 
tic time at scale r in the inertial range is 0(r/v r ) where 
v r = (r/rjY^Vjj and v v = (eiy) 1 ^ 4 . The eddy damping 
rate, which measures the rate of cascade to smaller r, is 
0{v r /r). Scales in the inertial range see only the eddy 
damping and do not see v. A physical cascade step must 
be more than a factor of 2 in 1/r; there is still substantial 
energy transfer at a factor of 4. The maximum number 
of effective cascade steps was discussed in Sec. 5 and will 
be discussed further in Sec. 7. 

The skewness (it 3 ) / (u 2 ) 3 / 2 is independent of r in the 
inertial range of K41. It measures the rate of energy cas- 
cade. There is no dependence on R\. Experiment and 
DNS indicate that the skewness is very nearly indepen- 
dent of r and R\ in actual intermittent turbulence. The 
same conclusion comes from several models of intermit- 
tency, including the beta model and multifractal models 
[15]. 

We noted in Sec. 3 that the energy cascade cannot be 
captured if the kinetic energy, which is a sum of squares 
in k or wavelet bases, is used in the Tsallis formalism. 
In (21), skewness is given to u by adding linear and tri- 
linear terms to the actual kinetic energy. This procedure 
requires examination. 

Suppose that 7 in (12) is non-stochastic but a is 
stochastic and constant in time. Then the PDF of y 
is obtained by a weighted integration over the variance 
parameter of a Gaussian PDF. Very wide classes of PDF 
can be generated in this way. See Sec. 8 for further dis- 
cussion. A Tsallis PDF with q > 1 is generated by % 2 
statistics for 1/cr 2 , as described in Sec. 4. No choice of 
weight function can give y a skewness. To obtain skew- 
ness, the weighted white noise o~w{t) must be replaced 
by a skewed forcing with finite correlation time. This is 
taken as the f(t) of (20). 

We have difficulty accepting the relation between tur- 
bulence and these manipulations with singlc-dcgrce-of- 
freedom models. Problems with characteristic times have 
already been noted. The inclusion of skewness in the 
models seems to be motivated by the empirical fact that 
the PDF of u is skewed, rather than any dynamical link 
to turbulence structure. 

In later work [16], (21) is modified and the dependence 
on R\ is dropped; cj\/R\ is replaced by c^/yf and it is 
stated that a match to experimentally measured skew- 
ness requires y/yf to be about 0.1. It is stated that c 
may vary with the experiment and may depend on q. 
Also, it is implied that the variation of q with r needs to 
be determined a posteriori. An additional parameter a 
is introduced yielding the effective energy 
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E{u) = ±\u\ 2a - c^F sign(u)(M Q - ±\u\ 3a ). (23) 

It is stated that a < 1 corresponds to a fractal phase 
space, in analogy to the beta-model picture of active ed- 
dies confined to a fraction of space. 

The complicated theoretical apparatus thus con- 
structed can give matches to the normalized PDFs of u 
measured in the central region of a Taylor-Couctte flow 
at macroscale Reynolds number R = 540 000. The flow is 
between a rotating inner cylinder and a stationary outer 
cylinder. The best fits are obtained by the empirical 
procedure of taking a = 2 — q, c^fyr = 0.124(<j — 1), 
and then choosing q to give minimum mean-square de- 
viation from the measured PDF. Once these expressions 
for a and c^/jf are fixed, the only free parameter is q. 
If the parameter a is not invoked, the range of l/(q — 1) 
needed for best fits to experimental data increases and 
there is some loss of accuracy in the best fits. The range 
of q — 1 used in fitting the measured curves is from 0.055 
at r/rj = 14400 to 0.168 at r/rj = 11.6, corresponding to 
a range of l/(q — 1) from 18.2 to 5.95. If the parameter 
a is not employed, the range increases. 

As already noted, the velocity difference u does not 
belong to a valid representation basis for the flow. Anal- 
ysis with given q and (3 must apply to appropriate modes 
in a valid underlying basis. The implication is that the 
given q and (3 must describe all statistics dominated by 
those modes, including but not limited to the statistics 
of transverse velocity differences. 

The modeling that yields (23) has been refined [17] 
by appeal to extended self-similarity (ESS) to determine 
(3 and thereby reduce the number of free parameters. 
Close matches to experimental estimates of scaling ex- 
ponents are obtained. An effort in a different direction 
[18] finds that the Tsallis entropies of subsystems can be 
made quasi-additive by adjustment of the q values for 
the subsystems. This leads to another formula for the 
variation of q with r. 

We have two principal difficulties with models in which 
q varies with r. The first is that the analysis appeals to 
fluctuations of effective equilibrium temperatures. Equi- 
librium temperatures, fluctuating or not, are inconsis- 
tent with the dominance, in the turbulence dynamics, of 
the strong and on-average unidirectional multistep en- 
ergy cascade. The second difficulty, noted in Sec. 3, is 
that modes whose scales differ by less than a decade have 
substantial dynamical coupling so that assigning them 
different q values seems unjustified. 

7. MULTIFRACTAL ANALYSIS OF VELOCITY 
FLUCTUATIONS 

Arimitsu and Arimitsu [3] use the relation of Tsallis 
statistics to multifractal processes to develop a multi- 
fractal model of high Reynolds number isotropic turbu- 



lence that gives detailed predictions of scale-dependent 
statistics. In addition, scaling exponents of velocity dif- 
ferences and other properties are derived and compared 
with the approximately isotropic DNS of [13]. Other, em- 
pirically developed and tuned, multifractal models have 
previously been found to agree well with a variety of mea- 
surements [19]. 

As noted in Sec. 5, the starting point of the analysis 
by Arimitsu and Arimitsu is the relation (15) between 
q and the extremes of the multifractal spectrum f(a) of 
a one-dimensional nonlinear map that displays sensitiv- 
ity to initial conditions. This relation is generalized and 
then applied to the multifractal spectrum of a turbulence 
model. An explicit value for q is found in terms of the 
intermittency exponent ji of the model. 

Some remarks on the physical nature of the inertial- 
range cascade: Viscosity is not felt directly at a scale 
r well within the inertial range. Instead the cascade of 
kinetic energy is characterised by a transfer rate n(r). 
n(r) is not positive-definite, and the DNS shows that 
it has a finite probability of negative values [13]. In 
contrast, all spatial averages of viscous dissipation arc 
positive. This contrast is recognized in [19]. Finite time- 
scales are associated with the cascade process at the scale 
r; there is not an instantaneous relation between cascade 
rate over a domain of size r and dissipation averaged over 
the domain. Moreover the energy cascade is sufficiently 
diffuse that there is significant transfer at scale ratios 4 
and larger. 

Arimitsu and Arimitsu use the multifractal formalism 
described in Sec. 5 to fit scaling exponents of velocity dif- 
ferences presented in [13]. Both longitudinal and trans- 
verse differences are fitted. They find that an optimum 
fit over a wide range of space intervals r is obtained by 
taking ji = 0.327 and n values that range from 21.5 to 
4.00 for r/rj values that range from 2.38 to 1220. 

Maximum admissible n values for various r/rj values 
follow from the discussion in Sec. 5. We noted there 
that the effective number of grid point spacings at the 
beginning of the unforced cascade range is ~ 209; that 
there is substantial scale nonlocality in the turbulence 
coupling; and that there must be a non-asymptotic re- 
gion at the beginning of the scaling range. Since rj in the 
DNS is approximately 1 grid point spacing, this leads to 
the conclusion that the maximum admissible n value for 
r w i] is less than 9. Extension of this argument to other 
r values gives n = 8, 7, 6, 5, 4, 3, 2, 1, as upper limits 
for r/rj = 1.95, 3.91, 7.81, 15.6, 31.3, 62.5, 125, 250, 500. 

All of the n values used in the optimum fitting in [3] 
greatly exceed the maximum admissible n values given 
above. We are thereby stymied in considering further the 
fits to velocity-difference statistics and scaling exponents 
of structure functions. We cannot assess the proposal 
that there are two multifractal scaling ranges separated 
by a crossover value of r. 
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8. CONCLUDING REMARKS 

The inviscid conservation of kinetic energy in incom- 
pressible 3D NS flow gives important information about 
the energy cascade process at large Reynolds number and 
the associated skewness of the velocity field. Energy con- 
servation fails to provide information about higher statis- 
tics - in particular, the development of intermittency in 
the cascade. A number of mechanisms that affect higher 
statistics are known. Vortex stretching implied by the 
Kelvin circulation theory and the existence of several 
kinds of instability play important roles. However, it is 
not known how to put all these things together into a de- 
ductive theory of turbulence statistics. Instead, there are 
a variety of models of higher statistics that have meager 
or nonexistent deductive support from the NS equations 
but can be made to give good fits to experimental mea- 
surements. 

A straightforward application of Tsallis statistics 
to turbulence suffers the same basic limitation as 
Boltzmann-Gibbs statistics: the equilibria are based 
solely on the kinetic energy of the flow. It was noted 
in Sec. 3 that all the degrees of freedom enter symmetri- 
cally into the kinetic energy and therefore are equally ex- 
cited in equilibrium, whether Boltzmann-Gibbs or Tsal- 
lis. Two strategies that have been used to break the 
symmetry are the construction of effective modal ener- 
gies that differ from the actual kinetic energy and the 
use of different values of the nonextcnsivity parameter 
q for different parts of the system. Some of the work 
that uses these strategies is described in Sees. 5-7. As 
explained in Sec. 3, the result is no cascade of energy 
between groups of modes at different scales, whether or 
not skewness appears in the PDFs. 

An absurd, but perhaps illuminating, manipulation of 
effective energies that breaks the symmetry over modes 
can be performed within Boltzmann-Gibbs statistics. 
Any desired PDF P n (y n ) can be realized for a mode y n by 
taking -j3e n (y n ) = \nP n (y n ) in (6) to determine e n (y n ). 
Pn(y n ) can be constructed to include skewness, but again 
the result does not describe the energy cascade. Energy 
transfer requires nonzero correlations of form {y n ymVr), 
where n,m,r label modes with differing wavenumbers or 
other scale labels. 

Beck and his coworkers treat velocity differences across 
distances r as dynamic flow variables, which they are not. 
A problem with their approach is that each velocity dif- 
ference is modeled by a modified Langevin equation with 
white-noise forcing. This conflicts with the fact that the 
driving of modes of given spatial scale is dominated by 
modes of larger scale that have finite and longer charac- 
teristic times. 

Arimitsu & Arimitsu generalize single-variable multi- 
fractal attractors for which Tsallis statistics is exact and 
propose a link to multifractal turbulence models. A key 



quantity in the fits to DNS statistics is the number of 
nominal "eddy size" halvings. Impressive fits are made if 
this number is freely disposable, but they are destroyed 
if the number is limited to meaningful values. 

The PDFs that have been fitted in applications of 
Tsallis statistics to NS turbulence have benign shapes: 
smooth and nearly symmetrical with flaring skirts. Such 
shapes can be generated in many ways. For example, 
consider the fits to acceleration PDF examined in Sec. 
5. With all questions of consistency and validity put 
aside, the fact is that both Beck and Arimitsu & Arim- 
itsu present fits to the entire range of the DNS PDF. Yet 
the fitting algorithms and physical models used by these 
authors are very different. Fits to PDFs of velocity incre- 
ments, and scaling exponents, have been obtained from 
a variety of models constructed by Beck. More challeng- 
ing would be shapes like the strongly asymmetrical and 
sharply cusped PDF of velocity gradient in Burgers tur- 
bulence forced at low wavenumbers [20] . 

Beck and Cohen [21] consider general ensembles of 
Boltzmann-Maxwell distributions over which is itself 
stochastically distributed. The j3 distributions explicitly 
examined include x 2 and log-normal. If the energy is a 
sum of squares of modal amplitudes, the PDFs of mode 
amplitudes over such ensembles can be given protean 
shapes, with a smooth central peak and flaring skirts, 
by adjusting the (3 distribution. A wide variety of empir- 
ical PDFs can be fitted in this way. In order to make an 
application to cascading turbulence, the interaction of a 
mode like, say, a k or wavelet amplitude with the rest 
of the system would have to be modeled as an interac- 
tion with a heat bath whose /? fluctuates. An ensemble 
of equilibria with heat baths is inconsistent with the en- 
ergy cascade, in which there is a strong, on-average uni- 
directional, flow of energy from larger to smaller scales. 
Strong here means that the cumulative mean energy flow 
through, say, a double octave of wavenumbers in one lo- 
cal characteristic time is the order of the mean energy in 
the double octave. 

Multifractal cascade models raise the general issue of 
distinction between what is descriptive of physical be- 
havior and what can be used for analysis of data. In 
[19] experimental data are analyzed to yield generalized 
dimensions D p and the corresponding multifractal spec- 
trum over the range — 20 < p < 20. Then the multifractal 
model is tuned to give an optimum (and excellent) fit to 
the data. Multifractal models may or may not express 
well the cascade physics at large but finite Reynolds num- 
bers. 
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